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A graphical  boundary  matching  procedure  was  used  to  determine  a likely  disk/blade  interface 
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disk  fillet,  relative  to  the  average  stress  in  the  neck  section  of  the  disk.  In  addition,  the  SCF 
was  found  to  decrease  with  decreasing  interface  friction.  Parametric  modification  of  the  disk 
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\ 20.  Abstract  (Continued) 

^ times  the  present  one.  For  all  studied  combinations  of  interface  force  distribution,  coefficient  of 
friction  and  fillet  geometry  the  point  of  maximum  tensile  stress  occurred  farther  into  the  fillet 
(away  from  the  contact  region)  than  the  point  at  which  cracks  appear  to  initiate  (near,  or  just 
into  the  inboard  edge  of  the  contact  area).  However,  different  combinations  of  interface  force 
and  fillet  geometry  result  in  different  degrees  of  spread  of  the  tensile  stress  concentration  from 
the  fillet  area  into  the  contact  area.  Finally,  the  effects  of  some  assumed  interface  force  per- 
turbations due  to  non-radial  (non-centrifugsJ)  loading  are  discussed. 
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NUMERICAL  PARAMETRIC  STRESS  ANALYSIS  OF  THE  TF-30 
TURBINE  ENGINE  THIRD-STAGE  FAN-BLADE/DISK 
DOVETAIL  REGION 


INTRODUCTION 


1 This  report  details  the  numerical  portion  of  a com- 
bined experimental /numerical  study  of  the  stresses  in  the 
third  stage  fan-blade/disk  assembly  of  the  TF-30  turbojet 
engine  used  in  the  Flit  aircraft  (see  Fig.  l),  A discussion 
of  the  background  for  this  study,  as  brought  on  by  crack- 
induced  in-service  failures  of  the  TF-30  engine,  can  be 
found  in  NRL  Report  8li*9,  which  contains  the  results  of  the 
experimental  program  [l].  The  experimental  study  of  both 
“he  disk  and  blades  involved  their  response  to  a single 
interface  loading  condition,  as  produced  by  an  axial 
(radial)  load  applied  to  the  blades  in  the  test  setup  of 
Fig.  2.  The  numerical  studies  reported  here  were  concen- 
trated on  the  disk/lug  contact/fillet  area,  and  involved 
the  sensitivity  of  this  region  of  the  lug  to  a variation 
of  interface  loading  parameters  such  as  contact  force 
distribution  and  normal-f orce / shear ( fri ct ion ) -for ce  ratios. 


2.  SIMULATION  SYSTEM  DESCRIPTION 


The  parametric  finite  element  analysis  was  performed 
using  a newly  developed  interactive,  graphics-oriented 
package  called  the  Two-dimensional  O_rthotropic  Real-Time 
Analysis  Library  (TOTAL)  [2].  It  consists  of  a set  of 
modular  programs  linked  by  a broadly  structured  disk-file 
data  base.  Each  of  the  thirty  programs  implements  a speci- 
fic task  in  the  simulation/solution/interpretation  proce- 
dure, and  creates  data  files  to  be  used  as  needed  by  other 
programs  further  along  the  solution  path.  All  of  the  pro- 
grams are  auto-prompting  and  to  a considerable  extent  auto- 
instructional  with  all  of  the  intricacies  of  data  file 
manipulation  hidden  from  the  user.  The  use  of  modular  pro- 
grams is  intended  to  encourage  numerical  experimentation 
(or  parametric  analysis).  For  example,  complete  or  partial 
changes  in  files  containing  boundary  conditions,  material 
properties  or  geometric  information  can  be  quickly  made  by 
programs  made  to  perform  specific  modification  tasks. 

Note:  Manuscript  submitted  December  1,  1977. 
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The  problem  can  then  be  resolved  without  reentering  those 
quantities  which  did  not  change.  Similarly,  compatible 
result  files  can  be  linearly  combined  to  study  the  effect  of 
superposed  solutions.  These  features  were  utilized  exten- 
sively in  this  study. 

3.  GENERAL  SOLUTION  PROCEDURE 

The  primary  mode  of  failure  for  the  third  stage  of  the 
TF-30  engine  has  been  observed  to  be  a fracture  in  the  disk- 
lug  contact/fillet  area,  with  some  similar  but  less  frequent 
fracture  of  the  blade.  Therefore,  the  numerical  simulation 
portion  of  the  TF-30  stress  analysis  concentrated  on  the  lug 
area  by  attempting  to  answer  the  following  three  questions. 

1.  What  are  the  most  likely,  and  most  extreme, 
contact  loading  conditions  expected  at  the  blade/lug 
interface? 

2.  What  are  the  maximum  values,  and  distribution 
characteristics,  of  the  tensile  stress  field  in  the  contact/ 
fillet  areas  of  the  present  disk-lug  configuration,  when 
subjected  to  the  above  loadings? 

3.  Are  there  minor  geometric  modifications  which 
will  decrease  the  tensile  stresses? 

The  program  of  studies  chosen  to  answer  these  questions 
consisted  of  the  following  steps. 

1.  Choice  of  suitable  general  geometry  and  boundary 
conditions  for  the  various  numerical  models  to  be  used, 
based  on  radial  (centrifugal)  loading  of  the  blade. 

2.  Preliminary  analysis  of  coarse-mesh  blade  and  lug 
models,  loaded  separately  and  also  in  a linked  (welded 
along  contact  surface)  condition. 

3.  Study  of  the  effect  of  distribution  of  interface 
forces,  and  of  the  relative  mix  of  normal  and  shear  (fric- 
tion) forces. 

U.  Choice  of  the  most  likely  distribution  of  inter- 
face forces  based  on  the  compatibility  of  interface  dis- 
placements for  various  assumed  force  distributions. 

5.  Study  of  the  effect  of  fillet  geometry  modifi- 
cation in  the  lug  using  a refined  mesh  and  a reference 
interface  loading. 
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ge:iehal  model  characteristics 


Basic  Geometries 


Due  to  the  relatively  small  angle  (5°)  between  blade 
and  lug  center  lines,  and  also  the  restriction  to  rectang- 
ular cartesian  constraints  with  the  TOTAL  system,  a rec- 
tangular (X-Y)  modeling  was  chosen,  rather  than  a cylindri- 
cal (R-0)  one.  The  region  to  be  modeled,  taking  advantage 
of  symmetry  along  the  center  lines  of  both  blade  and  lug, 
is  illustrated  in  Fig.  2.  The  choice  of  a rectangular 
modeling  had  the  added  advantage  of  providing  for  comparison 
with  the  results  of  the  experimental  photoelastic  analysis. 
The  setup  for  the  experimental  analysis  is  shown  in  Fig.  3. 
In  order  to  facilitate  comparison,  the  boundaries  of  the 
numerical  model  were  taken  directly  from  the  working  draw- 
ings of  the  photoelastic  model,  except  for  the  elimination 
of  excess  material  at  the  top  of  the  blade  and  bottom  of 
the  lug.  The  effect  of  shortening  the  base  of  the  lug  is 
shown  to  be  negligible  later. 

The  mesh  of  Fig.  Ua  is  relatively  coarse,  and  was  used 
primarily  for  the  solution  of  the  two  parts  in  a linked 
condition,  and  for  the  few  individual  solutions  run  for  the 
blade.  For  the  more  extensive  studies  of  the  lug,  the 
refined  mesh  of  Fig.  5a-  was  used. 

Reference  Boundary  Conditions 

Applied  load:  The  assumed  physical  loading  is  a verti- 
cal (radial)  force  applied  to  the  blade,  as  in  the  experi- 
mental setup  of  Fig.  3.  This  loading  is  assumed  to  produce 
identical  force  distributions  on  all  contact  surfaces. 

(The  effect  of  non-radial  loading  components  is  briefly 
discussed  in  the  Appendix).  For  the  reference  loading  at 
the  interface  of  each  part,  the  following  force  systems, 
each  uniformly  distributed,  were  chosen: 

1.  A uniformly  distributed  unit  normal  force,  applied 
at  the  seven  nodes  defining  the  interface. 

2.  A .175  (coefficient  of  friction)  distributed  shear 
force,  applied  at  the  same  seven  nodes.  (Choice  of  this 
friction  coefficient  is  based  on  initial  interpretation 
of  the  photoelastic  analysis,  and  on  previous  numerical 
analyses  of  similar  configurations.  Later  analysis  of 
experimental  data  indicated  a friction  coefficie-nt  of  .12  to 
.16). 

Note:  for  comparison  purposes  results  for  all  loading 
cases  and  geometries  are  normalized  with  respect  to  the 
combined  vertical  (Y)  component  of  this  reference  loading. 


3 


I 


Constraints:  For  the  lug,  the  left  and  right  vertical  edges 
were  fixed  in  the  X-direction  and  free  in  Y,  by  invoking 
free  in  symmetry  conditions.  The  bottom  edge  was  fixed  in 
Y and  X as  a first  approximation  to  the  condition  in  the 
real  part  and  in  the  photoelastic  model.  For  the  blade,  the 
vertical  edges  were  similarly  fixed  in  X,  and  the  top  edge 
was  constrained  in  Y. 

Linked  condition:  For  the  solution  of  the  two  parts 
linked  at  the  interface,  the  X-constrairt  on  the  vertical 
edges  and  the  Y-constraint  on  the  bottom  edge  is  maintained, 
and  a uniform  distributed  Y-force  is  applied  to  the  top  edge 
of  the  blade.  (See  Fig.  3). 

5.  PRELIMINARY  ANALYSIS  AND  MODEL  EVALUATION 


Note:  Although  the  quantity  of  principal  interest  in 
this  study  is  the  tensile  stress,  comparisons  of  various 
models  in  this  section  are  made  on  the  basis  of  isochro- 
matic  (maximum  shear)  patterns,  since  they  afford  an  imme- 
diate comparison  with  photoelastic  results.  In  addition 
some  isopachic  (contours  of  summed  normal  stresses,  a + o^ ) 
fringes  are  also  plotted  for  comparison  with  holographic 
patterns.  At  a free  boundary  both  of  these  patterns  reduce 
to  the  tangential  stress. 


The  results  of  this  section  are  best  presented  by  dis- 
cussing Figs,  h and  5.  In  Fig.  La-c  the  overall  maximum 
shear  pattern  is  shown  and  also  a local  enlargement  in  the 
critical  area  (with  fringe  numbers)  for  the  coarse  lug  model 
with  the  combined  reference  loading  and  constraints.  The 
resulting  stress  concentration  factor  (SCF)  is  U.8. 


If  the  same  reference  boundary  conditions  are  applied  to 
the  refined  lug  model.  Figs.  5b-c  show  that  the  general 
characteristics  of  the  shear  pattern  do  not  change,  but  that 
there  is  some  change  in  the  critical  area  (approximately 
1/2  fringe)  due  to  the  greater  resolution.  The  SCF  in- 
creases W to  5.2. 


Figs.,  6a-b  are  a similar  representation  of  what  happens 
if  the  base  of  the  unrefined  model,  with  the  reference 
loading,  is  extended  in  a manner  similar  to  the  photoelastic 
'model.  It  is  seen  that  there  is  little  effective  change,  and 
none  in  the  critical  region. 


r The  last  perturbation  on  the  basic  model  is  illustrated 

^ in  Figs.  7a-b.  The  reference  interface  loading  of  the 

refined  mesh  is  maintained,  but  the  right  vertical  edge  is 
now  fixed  in  both  X and  Y.  Again,  although  there  is  some 
considerable  change  of  the  shear  pattern  near  the  right 
edge,  the  critical  area  remains  unchanged. 
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For  comparison  purposes,  the  photoelastic  isochromatic 
and  isopachic  fringe  patterns  are  shown  in  Figs.  3 and  9* 

Also  the  numerical  isopachic  patterns  for  the  reference 
loading  of  the  coarse  and  fine  lug  models  are  given  in 
Figs.  10a  and  10b. 

As  a final  task  in  the  preliminary  analysis,  the  combined 
blade/lug  configuration  of  Fig.  2 was  loaded  as  described 
above,  and  the  properly  normalized  shear  contours  plotted  in 
Figs,  lla-b.  For  the  combined  geometry,  and  for  the  critical 
area  of  the  lug,  it  is  seen  that  this  physically  unrealizable 
case  has  global  features  similar  to  the  preceding  cases. 
However,  the  critical  areas  for  both  the  blade  and  the  lug 
are  so  close  to  the  interface  that  this  linked  condition 
cannot  be  expected  to  produce  a stress  concentration  similar 
to  that  of  the  two  parts  loaded  simultaneously  and  in  con- 
tact . 

The  overall  conclusions  to  be  drawn  from  these  prelim- 
inary investigations  are  that  the  refined  mesh  provides 
enough  improvement  over  the  coarse  one  to  Justify  the  small 
extra  computation  cost,  that  a deep  lug  base  is  not  required, 
and  that  the  vertical  edges  are  best  left  free  in  Y but 
fixed  in  X. 

6.  INTERFACE  FORCE  VARIATION 

In  order  to  study  the  effect  of  location  and  distri- 
bution of  the  interface  force-s  on  the  magnitude  of  the  stress 
concentration  in  the  lug,  a series  of  fourteen  distributed 
and  concentrated  loadings  were  applied  to  the  contact  sur- 
face of  the  refined  lug  model.  Seven  normal  loadings,  each 
with  a resultant  force  equal  to  that  of  the  reference  nor- 
mal loading,  are  shown  in  Fig.  12.  The  other  seven  loadings 
consisted  of  corresponding  shear  forces  with  the  same  distri- 
bution and  location  as  the  illustrated  normal  loadings.  The 
total  resultant  shear  force  in  each  case  was  equal  to  that 
of  the  reference  shear  loading,  or  equivalently  17.5^  of  the 
corresponding  normal  force.  It  should  be  mentioned  here 
that  single  point  forces  are  included  not  as  likely  physical 
situations,  but  simply  as  hypothetical  extreme  bounds  on  the 
parametric  analysis. 

For  ease  in  interpretation  of  results,  the  following 
procedure  has  been  used  in  naming  the  loading  cases.  The 
first  character  denotes  normal  (N),  shear  (S)  or  combined 
(NS)  loading.  The  second  character  denotes  distributed  (D) 
or  concentrated  (C).  The  third  character  denotes  location 
of  the  force  resultant  at,  or  near,  the  bottom  (B),  middle 
(M)  or  top  (T)  of  the  contact  surface  (or  combinations  of 
3,  M,  T) . 


The  results  of  this  parametric  analysis  are  presented  in 
the  table  of  Fig.  12,  which  contains  the  SCF  for  each  of  the 
fourteen  loads,  and  also  for  the  combination  of  correspond- 
ing normal  and  shear  loads.  In  each  case,  the  SCF  is  ob- 
tained by  dividing  the  maximum  fillet  tensile  stress  by  the 
average  tensile  stress  across  the  minimum  section  connecting 
the  fillet  area  with  the  right  hand  vertical  edge. 


Effect  o f d i s tributlon 


The  following  trends  are  evident  from  a study  of  Fig.  12. 
As  the  centroid  of  force  distribution  at  the  lug  interface 
shifts  from  the  bottom  towards  the  top  end,  the  normal  con- 
centration increases  but  the  shear  concentration  decreases. 
However,  for  the  assumed  reference  loading  ratio  (l.O  normal 
to  .175  shear)  the  combined  stress  concentration  increases. 
The  overall  concentration  is  bounded  by  a value  of  6.0  for 
a concentrated  load  at  the  top  of  the  interface,  as  com- 
pared to  5-2  for  the  reference  uniformly  distributed  loading 
(NDM  + SDM) . 


Effect  of  friction 


For  any  two  corresponding  normal  and  shear  loadings 
(e.g.,  NDB,  SDB),  the  effect  of  changing  the  friction 
coefficient  f,  while  maintaining  the  total  vertical  force 
constant,  can  be  found  by  substituting  values  from  the  table 
into  the  following  formula: 


SCF 


SCFS  + 


1 + ( .175  - f ) 


SCFN 


where  SCFS  and  SCFN  are  the  shear  and  normal  values  from  the 
table  of  Fig.  2.  As  an  example,  for  loadings  NDB  and  SDB, 
and  for  half  the  reference  friction  (f  = .0825),  we  get 
SCF  = U.8.  It  is  seen  from  the  form  of  the  above  formula, 
and  from  the  table  values,  that  a lower  coefficient  of 
friction  always  results  in  a lower  stress  concentration. 


7.  COMPATIBILITY  OF  INTERFACE  FORCES 


The  following  procedure  was  followed  in  estimating  the 
most  likely  interface  distribution  from  among  the  four 
distributed  load  cases  of  Fig.  12.  It  is  based  on  the 
assumption  that  the  true  interface  distribution  should  pro- 
duce, in  both  the  blade  and  the  lug,  at  least  the  following 
two  features : 


1.  The  same  rotation  of  the  interface  surfaces. 


2.  Maintenance  of  contact  along  the  interface  consis- 
tent with  the  assumed  interface  loading  distribution.  For 
example,  loadings  concentrated  midway  along  the 
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interface,  and  producing  a resulting  concavity  of  both  sur- 
faces, wouli  be  inconsistent. 

Since  the  nature  of  t. e chosen  interface  loading  var- 
iations (linear  variation  only)  could  not  be  expected  to 
produce  the  exact  interface  distribution,  a method  of  simple 
visual  comparison  of  the  approximate  compatibility  between 
displaced  interfaces  of  the  blade  and  lug  was  used.  In 
addition,  the  solutions  for  the  blade  were  obtained  using 
the  coarse  mesh  of  Fig.  2,  whereas  the  refined  model  of 
Fig.  5 was  used  for  the  lug.  It  was  not  felt  that  deter- 
mination of  the  blade  displacements  warranted  a refined 
modeling. 

As  an  initial  reference,  an  attem.pt  was  made  to  match 
the  contact  surfaces  of  the  distorted  geometries  of  the  two 
parts  for  the  reference  uniformly  distributed  loading.  It 
is  seen  from  Fig.  13a  that  the  difference  in  interface 
rotation  makes  it  impossible  to  match  the  two  end  points, 
and  that  there  is  an  inadmissable  concavity  of  both  surfaces. 

Next  the  assumption  was  made  that  contact  stresses  were 
concentrated  towards  the  lower  end  (loading  case  NDB  + SDS ) . 
The  distorted  geometries  resulting  from  this  loading  are 
shown  in  Fig.  13b.  It  is  seen  that  this  results  in  a much 
smaller  rotational  incompatibility,  but  that  the  concave 
curvature  is  still  present  along  both  surfaces. 


i Next,  contact  loading  concentrated  near  the  upper  end 

(NDT  + SDT)  as  used,  and  Fig.  13c  shows  the  resulting  over- 
\ layed  surfaces.  Hers  an  extreme  overlapping  at  the  opposite 

[ (lower)  end,  together  with  a concave  interface  curvature, 

is  observed. 

These  results  indicate  that  any  loading  spanning  the 
midpoint  of  the  contact  surface  produces  an  inadmissable 
concavity,  and  that  contact  stresses  should  be  distributed 
mors  towards  the  lug-fillet  end  of  the  interface.  This 
' suggests  the  use  of  a loading  which  is  zero  at  the  center 

S of  the  interface,  and  which  varies  linearly  towards  each 

end,  with  a greater  intensity  at  the  lug-fillet  end  as  in 
loading  case  NLIK  in  Fig.  12.  The  overlayed  displacements 
f are  shown  in  Fig.  13d  for  a 1.3  to  0.7  ratio  of  the  inten- 

1 sities  at  the  two  ends  of  the  contact  area.  It  is  seen 

( that  the  concavity  has  practically  disappeared,  while  the 

|t  rotation  of  the  two  surfaces  remains  approximately  equal, 

p The  conclusion  to  be  drawn  is  that  this  loading  represents 

j a reasonable  approximation  to  the  contact  distribution. 
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This  is  in  agreement  vit.t  most  contact;  solutions  in  which 
interface  forces  cluster  towards  the  ends  of  the  contact 
surface . 

8.  FILLET  MODIFICATIOn 

Since  the  small  fillet  radius  in  the  critical  region 
appeared  to  be  a prime  contributor  to  the  stress  concen- 
tration, a local  geometry  modification  was  undertaken  as 
follows.  The  existing  fillet  radius  (.055  in..)  was  system- 
atically increased  and  the  neighboring  fillet  at  the  bottom 
of  the  key-way  was  shifted  without  change  in  radius  in  order 
to  maintain  continuity  of  the  boundary.  Three  different 
fillet  radii  were  chosen.  These  three  corresponded  to 

multiples  of  2,  3 and  U times  the  origina^  radius.  This 
modification  is  best  understood  by  looking  at  Figs.  lU-lT. 

A fourth  modification  consisted  of  replacing  the  two 
fillets  and  their  connecting  vertical  boundary  with  the 
largest  possible  single  circular  arc  which  is  tangent  to 
both  the  ^5°  contact  boundary  and  the  horizontal  boundary 
at  the  bottom  of  the  keyway,  without  intruding  into  the 
contact  area.  This  geometry  is  shown  in  Fig.  IT,  and  is 
seen  to  be  quite  similar  to  Fig.  lU . 

The  reference  loading  was  applied  to  all  four  modified 
geometries.  The  resulting  tensile  stress  distribution  in 
their  respective  critical  areas  is  plotted  for  each  geo- 
metry in  Figs.  lU-lT.  These  figures  also  contain  their 
respective  SCF. 

It  is  seen  from  these  stress  distributions  that  increas- 
ing the  fillet  radius  results  in  a decrease  in  the  stress 
concentration  and  a shifting  of  the  critical  area  away  from 
the  contact  surface.  However,  the  improvement  reaches  its 
optimum  in  Fig.  15  (3XR),  and  the  stress  then  starts  to 
rise  with  increasing  radius.  The  conclusion  to  be  drawn 
from  this  analysis  is  that,  as  expected,  the  stress  concen- 
tration can  be  decreased  up  to  2T%  by  increasing  the  fillet 
radius  to  an  optimum  value. 

9.  CONCLUSIONS 

On  the  basis  of  a two  dimensional  parametric  finite 
element  simulation,  the  following  conclusions  can  be  drawn 
about  the  value  and  distribution  of  the  edge  tensile 
stresses  in  the  critical  contact/fillet  area  of  the  TF-30 
third-stage  turbine  disk. 

1.  A likely  interface  condition  consists  of  very 
little  force  at  the  center,  with  approximately  twice  as  much 
force  at  the  lug-fillet  (inboard)  end  as  at  the  other  end. 
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2.  T;ie  stress  20  nc  en  t r at  i on  factor  (2C?)  in  the  lu^r 
fillet  for  the  f.ost  likely  interface  distribiatior.  and  for  a 
friction  coefficient  of  .175  is  5.3. 

3.  The  worst  case  SCF  is  6.0. 

U.  Decreasing  the  coefficient  of  friction  decreases 

t h e S C ? . 

J 

5.  Increasing  the  fillet  radius  can  decrease  the 
SCF  by  up  to  27^. 

.6.  In  all  cases  the  point  of  maximum  SCF  occurs  f 

farther  into  the  fillet  (away  from  contact  area)  than  the 
point  at  which  cracks  have  been  observed  to  start. 

7.  The  location  of  maximum  SCF  and  general  distri- 
bution are  consistent  with  the  experimental  results  of  NRL 
Report  31U9. 

As  a final  graphical  conclusion  the  shear  and  isopachic 
stress  contours  for  the  most  likely  loading  are  shown  in 
Fig.  18,  and  linear  plots  of  the  tensile  stress  around  the 
fillet  and  along  uhe  direction  of  observed  fracture  are 
shown  in  Fig.  19.  It  should  be  noted,  from  comparison  of 
Figs.  18a  b and  Figs.  8 and  9,  that  the  "most  likely" 
interface  loading  does  produce  a shear  stress  pattern  which 
is  closer  to  the  experimental  pattern  than  any  other  distri- 
bution studied. 
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During  an  oral  presentation  of  the  hasic  results  on 
l!+  February  19T7,  discussions  with  XAVAIR  and  Pract  and 
Whitney  representatives  indicated  that  they  observed  crack- 
ing to  originate  fro";  a point  inside  the  blade / disk-lug 
contact  region,  or  removed  from  the  point  of  maximum  tensile 
edge  .stress  as  determined  numerically  and  experimentally. 
Subsequent  examination  of  the  full  disk/blade  assembly  pro- 
vided to  NRL  showed  the  following:  (l)  1^4  cracks  origin- 
ating from  Just  inside  the  bearing  area,  all  in  lugs  on  the 
same  side  of  their  respective  blades;  (2)  2 cracks  origin- 
ating Just  outside  the  contact  area  (just  into  the  fillet 
area),  both  in  lugs  on  the  side  of  the  blade  opposite  to 
the  other  fourteen.  This  type  of  consistent  unsynmetric 
behavior  seemed  to  indicate  a sensitivity  to  ncn-radial 
blade  forces  (assumed  smaller)  which  might  be  superimposed 
on  the  centrifugal  forces.  This  situation  led  to  the 
following  brief  consideration  of  the  effect  of  non-radial 
forces  as  they  might  cause  a consistent  lack  of  symmetry 
in  the  tensile  edge  stresses,  relative  to  the  discrepancy 
between  locations  of  maximum  stress  and  crack  origination. 
Since  the  approach  used  in  choosing  possible  perturbations 
in  interface  loading  caused  by  non-radial  blade  loading 
involved  a considerable  amount  of  conjecture,  the  results 
should  be  interpreted  in  this  context. 

The  initial  analysis  was  based  on  an  assumed  centri- 
fugal (radial)  loading  due  to  rotation  of  the  disk.  However 
there  is  also  a wind  force  component  perpendicular  to  the 
blade  axis  which  tends  to  cock  the  blade  to  one  side.  It 
is  evident  that  this  cocking  action  will  shift  the  contact 
forces  towards  opposite  ends  of  the  intended  contact  sur- 
face on  either  side  of  the  blade.  On  one  surface  this 
cocking  will  tend  to  drive  the  distribution  from  the  "most 
likely"  situation  (loading  case  NLIK  + SLIK,  Fig.  12) 
towards  a more  balanced  distribution  (loading  case  HD3T  + 
SD3T).  This  increases  the  SCF  less  than  2%,  but  does  not 
appreciably  effect  its  location.  On  the  other  surface  the 
contact  forces  will  be  driven  towards  the  lug  fillet, 
approaching  a concentrated  force  condition  (loading  case 
(!IC3  + SC3).  Rather  than  using  a concentrated  force  to 
represent  this  situation,  three  new  loadings  were  applied 
to  the  lug,  each  with  normal  and  shear  (17.5^)  resul- 
tants equal  to  the  reference  loading.  (The  "most  likely" 
loading  is  given  as  a reference  in  Fig.  20a,  and  the 
resulting  tensile  stresses  in  Fig.  20e.) 
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The  first  involved  a linear  variation  from  zero  at  the  center 
of  the  contact  area  to  a maximum  at  the  end  nearest  the 
fillet,  and  no  load  on  the  other  half  of  the  interface 
(Fig.  20h).  The  second  consisted  of  applying  the  full  load 
at  the  two  nodes  nearest  the  fillet,  with  the  closest  node 
getting  twice  as  much  force  as  the  next  one.  (Fig.  20c). 

The  third  consisted  of  shifting  the  two  loading  points  of 
case  two,  each  one  node  away  from  the  fillet,  leaving  the 
closest  node  unloaded.  (Fig.  20d).  This  last  loading 
case  is  intended  to  study  the  effect  of  a small  shift  of  the 
contact  area  away  from  the  fillet  area.  Plots  of  the  ten- 
sile stresses  in  the  f illet/contact  region  for  these  three 
loading  conditions  are  shown  in  Figs.  20f,  20g,  and  20h. 

It  is  seen  that  there  is  no  significant  change  in  the  value 
or  location  of  the  maximum  tensile  stress  in  relation  to 
the  reference  loading  case.  However,  Fig.  20g  shows  that 
shifting  the  contact  area  towards  the  fillet  induces  higher 
tensile  stresses  into  the  lower  region  of  the  contact  area. 
This  would  result  in  a greater  likelihood  of  cracks  pro- 
gressing from  the  lower  end  of  the  contact  surface  due  to 
the  presence  of  local  imperfections.  This  is  in  accord 
with  observations  of  a cracked  disc  in  which  all  cracks 
which  did  not  initiate  in  the  fillet  area  occurred  on  the 
contact  surface  for  which  the  blade-cocking  tends  to  shift 
the  contact  towards  the  fillet  region. 
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(b)  Dovetail  region 

Fig.  1 — TF-30  third-stage  disk  blade  assembly 
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(b)  Refined  mesh 


Fig.  10  — Numerical 
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(c)  Loading  NDT  + SDT  (d)  Most  likely  loading  NLIK  + SLIK 


Fig.  13  — Matching  of  interface  displacements 
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